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We perform a determination of the strong coupling constant using the latest ATLAS inclusive 
jet cross section data, from proton-proton collisions at ^/s = 7TeV, and their full information 
on the bin-to-bin correlations. Several procedures for combining the statistical information from 
the different data inputs are studied and compared. The theoretical prediction is obtained using 
NLO QCD, and it also includes non-perturbative corrections. Our determination uses inputs with 
transverse momenta between 45 and 600 GeV, the running of the strong coupling being also tested 
in this range. Good agreement is observed when comparing our result with the world average at 
the Z-boson scale, as well as with the most recent results from the Tevatron. 



I. INTRODUCTION 

The strong coupling constant as is one of the fun- 
damental parameters of the Standard Model of particle 
physics (SM). Testing the energy dependence (running) 
of as over a wide range provides an implicit test of QCD 
and probes potential effects from "New Physics" . 

In this paper we present a determination of as ex- 
ploiting the recently published ATLAS inclusive jet cross 
section data [1] , allowing for the first test of the running 
of as up to the TeV scale. The only other input based 
on experimental measurements consists of the NLO QCD 
proton parton density functions (PDFs), as obtained by 
various groups from other independent data sets. These 
PDFs allow us to compute, for a given ag value, the per- 
turbative NLO QCD prediction of the inclusive jet cross 
section, which is then corrected for non-perturbative ef- 
fects. Our evaluation procedure goes through the deter- 
mination of as(Af^) using the measured cross section in 
each data bin, followed by the propagation of the experi- 
mental uncertainties and bin-to-bin correlations. Finally, 
an unbiased averaging procedure is used and the uncer- 
tainties are propagated to the corresponding result. 

In Sec. II we discuss the information provided by the 
ATLAS inclusive jet cross section data, while in Sec. Ill 
we describe the theoretical prediction used in the as de- 
termination. Several possible procedures of data treat- 
ment are discussed in Sec. IV, while the corresponding 
results are presented in Sec. V. Our conclusions from this 
study are discussed in Sec. VI. 



II. INPUT DATA 

In this study we exploit the unfolded double- 
differential ATLAS inclusive jet cross section data, in 
bins of the absolute rapidity (|j/|) and as a function of 
the jet transverse momentum (px), together with the 
statistical and systematic uncertainties, and their cor- 
relations [1]. Jets are defined using the anti-fct jet algo- 



rithm [2] implemented in the FAST JET [3] package. This 
measurement, performed for anti-fct jets with a distance 
parameter R=0.4 and R=0.6 respectively, covers a large 
phase-space region going up to 4.4 in absolute rapidity 
and from 20 GeV to more than 1 TeV in p^ . The detailed 
analysis of the uncertainties and their correlations is one 
of the important achievements of this measurement. 

The dominant uncertainty in this measurement is 
due to the jet energy scale (JES) calibration. Other 
smaller systematic uncertainties are due to the luminosity 
measurement, multiple proton-proton interactions ("pile- 
up"), trigger efficiency, jet reconstruction and identifica- 
tion, jet energy resolution and deconvolution of detector 
effects (unfolding). The information on the (anti-) cor- 
relations between the various \y\ and pt bins is provided 
for these components, using a set of 87 independent nui- 
sance parameters. These nuisance parameters also pro- 
vide the information on the asymmetries of these uncer- 
tainties (mainly due to the steeply falling px spectrum) . 
This type of information cannot be incorporated in a co- 
variance matrix. Depending on the rapidity bin, the sys- 
tematic uncertainties account for 20 — 60% at lowpx, and 
20 — 40% at high px: while they are smaller for interme- 
diate jet transverse momentum values. 

The statistical uncertainties and their (anti-) correla- 
tions are provided through a covariance matrix for each 
rapidity bin. Here the correlations are due to the trans- 
fers of jets between the different px bins (performed in 
the unfolding), as well as to the fact that several jets 
going into the spectrum can be produced in the same 
event. The statistical uncertainties are in general much 
smaller than the systematics, except for the high px I'e- 
gion, where they can be larger. 

It has been pointed out in [1] that, when performing 
the comparison between data and the theoretical predic- 
tion, somewhat larger differences are observed for R=0.4 
than for R=0.6 jets. However, the x^/dof , computed tak- 
ing into account the correlations, is reasonable in both 
cases. In our study, we do not perform a simultaneous 
treatment of the two jet sizes, because the correlations 
of the uncertainties of these measurements (expected to 



be large) are not available. We use the measurement of 
jets with R=0.6 for determining our nominal result, while 
the comparison with the result using jets with R=0.4 pro- 
vides a systematic uncertainty. 



III. THEORETICAL PREDICTION 

Theoretical predictions for the inclusive jet cross sec- 
tions are calculated in pcrturbative QCD at next-to- 
leading order (NLO) using the NLOJET-f -I- [4] program, 
while the effects of hadronisation and underlying event 
activity on jet production have been studied using Monte 
Carlo generators. Each bin of the pcrturbative cross sec- 
tion is multiplied by the non-perturbative corrections to 
obtain the prediction at the particle level that can be 
compared to the data. 



A. Perturbative QCD calculations 

The inclusive jet cross section is calculated using the 
CTIO [5] NLO parton density functions. The pt and 
rapidity binnings follow the one from ATLAS inclusive jet 
cross section measurement [1] . For the sake of simplicity, 
the renormalisation /ip and the factorisation /ip scales for 
the central result are set to be the same in the calculation 

M = MR = /"F = pT""^ 

where p^'^'^ is the transverse momentum of the leading 
jet in a given rapidity bin. 

To allow for fast and flexible evaluation of the cross sec- 
tion using different values of as, as well as for evaluation 
of the PDF and scale uncertainties, the APPLGRID [6] 
software was interfaced with NLOJET-f + in order to cal- 
culate the perturbative coefficients once and store them 
in a look-up table. 



B. PDF inputs 

There are several PDF fitter groups who provide alter- 
native PDF sets using slightly different theoretical ap- 
proaches and various experimental data from different 
scattering processes as input. Usually PDF groups pro- 
vide the central (best-fit, given an as{M'^) value) PDF 
set together with a collection of uncertainty PDF sets. 
The central set is used for the theory prediction, while 
the uncertainty sets are needed for the evaluation of the 
PDF uncertainties in the calculation. 

The best-fit CTIO PDF set is obtained with as{Ml) = 
0.118. In order to facilitate the studies of the dependence 
of observables upon the strong couphng CTIO PDF is 
provided with the special list of PDF sets, obtained by 
scanning the as(Mf). The PDF sets in this list were 
refitted while the as(Af^) was varied in the [0.113, 0.123] 
range in steps of 0.001. 



TABLE I: The range of cis-scan and the best-fit as(M|) 
values for PDF sets used in the analysis. 



PDF set 


asiMlf^^' 


as(M|)™ 


as(M|)-- 


step 


CTIO 


0.118 


0.113 


0.123 


0.001 


MSTW 2008 


0.1202 


0.110 


0.130 


0.001 


HERAPDF 1.5 


0.1176 


0.114 


0.122 


0.001 


NNPDF 2.1 


0.1190 


0.110 


0.130 


0.001 



Given the as-scan, we have calculated the inclusive 
jet cross section (T"^''(as,pT, \y\) for each as-point in the 
scan. Such a function provides a one-to-one correspon- 
dence between the agiM"^) and the jet cross section val- 
ues in each (px, \y\) bin. For the values of as outside 
the scan range (between scan points) a linear extrapola- 
tion (interpolation) is used. 

In contrast to CTIO, the central MSTW 2008 fit [7] 
suggests that as(M|) = 0.1202 and MSTW group pro- 
vides the as-scan in the range [0.110,0.130]. For the 
NNPDF 2.1 [8, 9] the central PDF set is fitted with 
as(M|) fixed to 0.1190 and the scan is performed for 
as values from 0.110 to 0.130 . Finally, for the HERA- 
PDF 1.5 [10] an as scan has been performed, from 0.114 
to 0.122, with the best-fit value at as(M|) = 0.1176. All 
the PDF sets provide the as-scan with the same stepping. 
The nominal as(A'/|) values for the various PDF sets, as 
well as the ranges of the scans, are summarized in table I. 



C. Non-perturbative corrections 

The fixed-order NLO QCD calculations predict hard 
parton-level cross sections, which need to be corrected 
for non-perturbative effects when compared to data. Af- 
ter the hard scattering phase of pp-coUision partons emit 
radiation and evolve down to lower scales. At a very low 
scale, when strong coupling becomes very large, the non- 
perturbative phase of collision starts, where coloured par- 
tons form colourless real hadrons. This process is called 
hadronisation. There is also underlying event (UE) ac- 
tivity, which happens in parallel to the hard scattering. 
The UE consists of interactions of beam remnants as well 
as from multiple parton interactions. 

The hadronisation process leads to a decrease of the 
jet pt, while the UE adds extra particles to the jet 
and results in an increase of jet transverse momentum. 
The effects of both hadronisation and UE are most pro- 
nounced at low pt- Since the effect of UE on the jet 
Pt increases with increasing size of the jet (jet radius R) 
while the effect of the hadronisation process decreases, 
if the jet radius grows, two separate sets of corrections 
and uncertainties were derived for jets with R — 0.4 and 
R = 0.6. The R dependence of the contributions of non- 
perturbative effects to the jet transverse momentum was 
calculated in [11, 12]. It has been shown that the hadro- 
nisation contribution is proportional to —1/R, while the 



TABLE II: Set of tunes used for the different generators. 



Generator 



Tune 



PYTHIA 6.425 



PYTHIA 8.150 
HERWIG++ 2.5.1 



AUET2B LO** 

AUET2 LO" AMBTl 

AMBT2B CTEQ6L1 

Perugia 2010 

Perugia 2011 

4C 

UE7000-2 



IV. as EVALUATION PROCEDURE 

In this section we describe several possible procedures 
for determining as{M'^) from the inclusive jet cross sec- 
tion data. First we explain how we exploit the experi- 
mental information in individual (px; \y\) bins, together 
with the propagation of the uncertainties and bin-to-bin 
correlations. Then we discuss several possible averaging 
procedures, for combining the statistical information of 
the various (px; \y\) bins, in order to obtain an average 
as{M^) value. We also discuss the specific properties of 
these various procedures. 



UE contribution varies as R^ . The leading R dependence 
was found to be quite similar for all jet algorithms. For 
R = 0.4, the correction factors are dominated by the ef- 
fect of hadronisation and are approximately 0.95 at jet 
Pt = 20 GeV, increasing closer to unity at higher px- For 
R = 0.6 the correction factors are dominated by the UE 
and are approximately 1.6 at jet px — 20 GeV decreasing 
to between 1.0 — 1.1 for jets above px = 100 GeV. 

The non-perturbative corrections are evaluated using 
leading-logarithmic parton shower generators. The cor- 
rections are derived using PYTHIA 6.425 with the AT- 
LAS AUET2B CTEQ6L1 tune [13]. These are evalu- 
ated as the bin-wise ratio of cross section calculated af- 
ter hadronisation and with simulation of the underlying 
event activity, to the one at parton-level and without un- 
derlying event. The parton-level cross section in each bin 
is then multiplied by the corresponding correction: 



A. Determination of as in individual (pt; \y\) bins 

An as value is obtained in each bin of the measure- 
ment, using the theoretical prediction and the experimen- 
tal cross section. The theoretical prediction provides, in 
each (px; \y\) bin, a one-to-one mapping between as and 
the inclusive jet cross section (see previous section). The 
nominal value of the measured cross section is mapped 
by this procedure to the nominal as value in the corre- 
sponding bin. All the experimental uncertainties of cross 
sections, together with their bin-to-bin correlations, are 
propagated to the determined as values, using pseudo- 
experiments (toys), as explained in the following subsec- 
tion. 



j2 particle 

dpT^d\y\ 



^2^NLO QCD 

dpT*d\y\ 



X <5"°"-P^'-* (p^^Myl) . (1) 



The uncertainty of the non-perturbative corrections 
is estimated as the maximum spread of the correc- 
tion factors obtained from PYTHIA 6.425 using the 
AUET2B LO", AUET2 LO**, AMBT2B CTEQ6L1, 
AMBTl, Perugia 2010, and Perugia 2011 tunes (PY- 
TUNE_350), and the PYTHIA 8.150 tune 4C [13- 
16], as well as the corrections obtained from 
HERWIG-F-h 2.5.1 [17] using the UE7000-2 [13] 
tune (see table II). In our analysis we have used the 
non-perturbative corrections evaluated in [1]. 

We would like to mention that an analytic estimate for 
the non-perturbative correction to the inclusive jet spec- 
trum was proposed [11]. It was used in [18] to compare 
the perturbative QCD calculations to the first ATLAS 
inclusive jet measurement [19]. The values for analytic 
and MC-based corrections were found to be compatible 
[18]. It would be very interesting to include the analytic 
evaluation of non-perturbative effects in the as deter- 
mination, once the flavour decomposition of inclusive jet 
spectra and UE parameters are measured experimentally. 



Propagation of data uncertainties and 
correlations 



A series of pseudo-experiments, each of them includ- 
ing a complete set of cross section values in all the (px! 
\y\) bins, are generated using the experimental input de- 
scribed in section II. 

The statistical fluctuations are generated using the in- 
formation provided in the published covariance matri- 
ces. Each matrix is diagonalised, and a set of inde- 
pendent random numbers (Gaussian distributions with 
widths given by the eigenvalues) is generated in the cor- 
responding space. The unitary matrices of eigenvectors 
are then used to propagate this set of random values to 
the physical space of the (px; \y\) bins. 

Each nuisance parameter for a given component of the 
systematic uncertainties is treated as 100% correlated be- 
tween all the bins. The various nuisance parameters are 
treated as independent with respect to each other. The 
distributions of the experimental uncertainties are mod- 
elled using asymmetric Gaussian shapes, centered at zero, 
with 50% probability of generating pseudo- measurements 
at negative or positive values. The two widths of the 
asymmetric Gaussian are equal to the corresponding un- 



certainties -'^. In a given pseudo-experiment, the contribu- 
tions of all the fluctuated nuisance parameters are added 
linearly to the nominal measured value. After adding 
the contribution of the statistical uncertainty, the value 
of the pseudo-measurement is obtained. 

For each pseudo-experiment, the same procedure as 
for the nominal measurements is applied to obtain an a^ 
value in each bin. This allows the propagation of all the 
experimental uncertainties, with their correlations and 
their asymmetric distributions, from the cross section 
level to the as values. 

In view of the study of various possible averaging pro- 
cedures, it is interesting to determine an approximate (af- 
ter the symmetrisation of uncertainties) covariance ma- 
trix of the as values. The effect of the needed symmetri- 
sation is rather small at this level, due to a partial com- 
pensation between the asymmetries of the cross section 
uncertainties and the non-linear effects in the theoretical 
prediction (for cr(as)). It is also interesting to determine 
a covariance matrix of the cross section measurements; 
however, the effect of the symmetrisation is larger at that 
level. 



C. Averaging procedures 

In order to study and avoid potential biases introduced 
by some averaging procedures, several options for per- 
forming the combination of the inputs in the various 
(pt; \y\) bins have been tested here. All of them aim 
at obtaining an average value of the strong coupling con- 
stant {af'^). 



the sense defined above. They also have the advantage 
that more precise contributions get larger weights. A 
variation of this procedure has also been tested, where 
the weights are proportional to the inverse of the squared 
statistical uncertainty. Actually, the statistical uncer- 
tainties exhibit much smaller correlations comparing to 
the systematic ones, motivating this additional test. 



3. Standard x^ minimisation 

The 3'"'^ averaging procedure consists of the minimisa- 
tion of a "standard" x^ with correlations 



X 



{a 



Nom 



a 



Av 



).c-i.( 



Nom 



a 



Av\T 



(2) 



Here, a^ ™ is the vector of nominal 



with respect to af'^ 
as values (obtained in individual (px; \y\) bins) entering 
the average, C is their covariance matrix, while a^ is a 
vector containing values equal to af"^ (to be fitted) and 
having the same size as a^"™. 

It has been shown in [20] that this approach, although 
very commonly used, can yield biased average values. Ac- 
tually, in presence of (not very well understood) strong 
correlations among the inputs, the average value can be 
even outside the range of the input values, which seems 
unacceptable. Indeed, the minimisation of the function 
in Eq. 2, with respect to a^ (equivalent to the minimisa- 
tion of the variance of the output of a weighted average, 
cf. Gauss-Markov theorem), yields (see for example [21]) 



i-c- 



(«^ 



1-C- 



JT 



(3) 



1. Simple average 

The first procedure consists of a simple average, where 
all the input as values have the same weight (equal to 
the inverse of the number of inputs). Although the final 
uncertainty on this average might not be optimal, the 
weights of the input as values are well behaved (i.e. they 
are all in the [0; 1] interval, with the sum equal to unity). 



2. Weighted average 

The 2"''^ procedure which has been tested consists of 
a weighted average, where the weights of the input as 
values are proportional to the inverse of their squared 
total uncertainties (obtained after symmetrisation). The 
weights used in this procedure are also well behaved, in 



^ This model of the distributions is motivated by the origin of the 
asymmetric uncertainties, due to the shape of the px spectra. 
However, as explained later, other possible models of these dis- 
tributions have been tested. 



where 1 is a vector with all the entries equal to unity 
and the same size as ag"". Just as in the two previous 
methods, the average resulting from the x^ minimisation 
is a weighted mean of the input values. However, the 
only constraint on these weights is that their sum equals 
unity. There is no constraint on the range of the weights 
and, in presence of strong correlations, these weights can 
even be negative or larger than unity. 

The explanation of the bias in this procedure, as pro- 
vided in [20] , makes responsible the input covariance ma- 
trix because of "the linearisation on which the usual er- 
ror propagation relies". However, two different problems 
of this approach have been pointed out in [22]. Actu- 
ally, Eq. 2 involves the inversion of the covariance ma- 
trix, which is ill-behaved (a small change in the covari- 
ance matrix can have an important impact on its inverse, 
hence on the weights and the average). Without special 
techniques, neither the inverse can be computed, nor the 
minimisation can be performed. Furthermore, this x^ 
definition treats all the uncertainties as absolute (i.e. cor- 
responding to additive effects), whereas at least some of 
them (like the uncertainty on the luminosity etc.) should 
be treated as relative uncertainties (i.e. corresponding to 
multiplicative effects). 



4- Modified x^ minimisation 



combination procedure uses the x^ test defined as 



A possible solution to these problems of the "standard" 
X^ minimisation has been provided in [20] . It consists in 
introducing, for each correlated systematic uncertainty, 
one nuisance parameter in the x^ definition, together 
with one extra constraint term. For a normalisation un- 
certainty, a scaling factor is applied to both data and the 
uncorrelated uncertainties. This is equivalent to intro- 
ducing the corresponding nuisance parameter as a scal- 
ing factor on the theoretical prediction. This alternative 
approach is used for example in [22, 23]. For an absolute 
uncertainty, the nuisance parameter acts a coherent ad- 
ditive shift of the theoretical prediction. However, it has 
been shown (see [24, 25]) that, if all the systematic un- 
certainties correspond to (are treated as) additive effects, 
the approach using nuisance parameters is equivalent to 
the one of the "standard" x^ with correlations (Eq. 2). 

This solution, treating the systematic uncertainties as 
nuisance parameters in the fit, is commonly known as 
"profiling". It often leads to reduced systematic uncer- 
tainties, as they are constraint by the x^ minimisation. 
However, this reduction strongly relies on the bin-to- 
bin correlations of the systematic uncertainties, as each 
nuisance parameter is treated as 100% correlated be- 
tween the bins. This strong assumption (exploited for 
propagating statistical information between the different 
phase-space regions) relies on information that is diffi- 
cult to check in a particular measurement and is there- 
fore, in most cases, not justified. This yields too op- 
timistic results, as the uncertainties on the shape and 
correlations for each component of the systematic uncer- 
tainties should also be taken into account. It should be 
pointed out that the same assumptions and uncertainty 
reductions are implicit in the "standard" x^ with correla- 
tions (Eq. 2). In order to better understand this problem, 
one can consider an example where, in the performance 
studies, a given component of the systematic uncertainty 
has been conservatively overestimated in some particular 
phase-space region. This conservative approach would 
be "exploited" by the fitting procedure, in order to re- 
duce the corresponding component in all the phase-space 
regions. 

In our case, the only systematic uncertainty for which 
the correlations are perfectly known is the one on the 
normalisation (i.e. the experimental uncertainty of the 
ATLAS luminosity). This motivates the fourth combina- 
tion procedure that we have tested, where the luminosity 
uncertainty is treated as nuisance parameter in the fit (it 
is "profiled"), while all the other uncertainties are in- 
tegrated in a partial covariance matrix. Therefore, the 
X^ definition will include a contribution from the com- 
parison between data and theory, involving this partial 
covariance matrix, plus a constraint on the luminosity 
shifts. 

As it is straightforward to propagate a luminosity vari- 
ation to the cross section, the first implementation of this 
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<^S ("s ) ■ PLJ 



Here, (ig " is the vector of nominal experimental cross 
section values, tTj^(a^^) is the vector of theoretical cross 
section values computed for a^^ (to be fitted), C_l;o- is 
the covariance matrix of the experimental cross section 
measurements, including all the uncertainties but the one 
for the luminosity, /3l is the luminosity scaling factor (ra- 
tio of the scaled to the nominal luminosity) , and cl is the 
luminosity relative uncertainty. The first contribution to 
the total x^ is given by a term (xl(/5l)) constraining 
the shifts of the luminosity nuisance parameter, while 
the second (x^ (cr (as); /3l)) comes from the comparison 
between data and theory at the cross section level. 

As explained in subsection IV B, the determination of 
the covariance matrix at the cross section level involves 
a larger symmetrisation, comparing to the determination 
of the covariance matrix of the as values. This motivates 
a second implementation of this combination procedure, 
driven by the x^ test: 



Xtot("s;/3L) 



/3l-1 



(5) 



{a 



Nom 



«r(/5L))-Cli^„,-(af"-af (/?l)) 



It involves the covariance matrix of the as val- 
ues (C_L;asj including all the uncertainties except the 
one for the luminosity), which is better defined, as ex- 
plained above, a^^ (/3l) is a vector containing values 
equal to a^^ (/3l) (the same for all the (px; |y|) bins) 
and having the same size as Q;g°"\ Here we have again 
the first contribution given by the term constraining the 
shifts of the luminosity nuisance parameter, while the 
second (xd ("s; /?l)) comes from the comparison between 
data and theory, at the a^ level. However, as the theo- 
retical predictions of the cross section values in the vari- 
ous (pt; |y|) bins exhibit slightly different Us dependen- 
cies (due to NLO QCD and non-perturbative effects) , the 
analytical dependence of a;^^ on a luminosity shift can 
be only approximate. In this second implementation, we 
use the following ansatz: 



Av 



(/3l 



1 



/3l 



1 



1.8 



(6) 



which is motivated by the small size of the non-linear ef- 
fects in the theoretical prediction, together with an em- 
pirical observation, as explained in section VE below. 



5. Geometrical average 

The first three combination procedures described 
above, make use of weighted averages. They distinguish 



between each other by the way how the weights of the 
inputs in the average are obtained. Another (the fifth) 
possibihty that we have considered for performing the 
combination consists in a geometrical average: 



Av 



n 



a 



Noni,i 



1 Nom.i 

where as 



runs over the n inputs of the combination. 



EXPERIMENTAL RESULTS 



In this section we present the results for the as deter- 
minations in the individual (px; \y\) bins, using anti-fct 
jets with R=0.6, followed by the results of the various 
combination procedures. We also discuss the issue of the 
disagreement between the results obtained for anti-fct jets 
with R=0.4 and R=0.6, respectively. 



A. Individual {pt; \y\) bins and their correlations 

Following the prescription described in section IV A, a 
nominal as value has been determined in each (px; |y|) 
bin, while the statistical and systematic uncertainties, 
together with their correlations, have been propagated 
using pseudo-experiments. Fig. 1 shows a few examples 
of as{M'^) distributions, as obtained from 10^ pseudo- 
experiments, in the central \y\ region, in four different 
Pt bins. Some asymmetries can be seen in these distri- 
butions. They are due to the asymmetries in the data 
uncertainties and to the non-linear dependence of the 
theoretical cross sections on as (from the theoretical cal- 
culation and its extrapolation to low and large as values, 
outside the region where the as scan is available). 

For very low and very large pt values, the total data 
uncertainties are relatively large. This is mainly due to 
the JES uncertainty, while the statistical uncertainties 
are also large at high px- In these regions of the px 
spectra, a fluctuation of the cross section by ±lcr (or 
more) yields as values for which the theoretical predic- 
tion is not reliable. Actually, the PDF determinations 
are not reliable for the corresponding large or low as 
values, while the perturbative calculations also exhibit 
problems at large as- It is for this reason that we dis- 
card the px bins with large experimental uncertainties in 
the final as average computation, and we estimate a sys- 
tematic uncertainty due to this choice. Because of their 
large uncertainties, these bins would in any case have a 
small impact on the average, therefore the estimated un- 
certainty associated to dropping these bins is small. A 
number of 42 (px; \y\) bins are used in the computation 
of the various averages. The corresponding ([px'";Pt'''^]; 
[Iz/lmin; |y|max]) raugcs, as wcll as the labels/numbers used 
for these bins, are given in table III. 



TABLE III: PT ranges i[p'T";pTl) used for the as deter- 
mination in various \y\ intervals ([|2/|min; |2/|max]). The third 
column indicates the label used for the \y\ bins, while the 
last column gives the range of numbers assigned to the corre- 
sponding (pt; 1 2/1) bins (ordered in increasing px and \y\). 



(7) [|2/|min;|y|max] [P' 



mm 
'T : 



max 

Pt 



](GeV) \y\ bin label bin numbers 



[0;0.3] 


[45;600] 


\y\i 


1-10 


[0.3; 0.8] 


[45; 600] 


\y\2 


11-20 


[0.8; 1.2] 


[45;500] 


|y|3 


21-29 


[1.2;2.1] 


[45;310] 


\vU 


30-36 


[2.1; 2.8] 


[60;210] 


\y\B 


37-40 


[2.8; 3.6] 


[80;110] 


|y|6 


41 


[3.6; 4.4] 


[80;110] 


\y\r 


42 



The same sets of pseudo-experiments have been used 
to determine the (total) covariance and correlation matri- 
ces, for the cross sections and as values respectively (see 
Fig. 2), in the (px; \y\) bins used for the final combi- 
nation. Although, as explained above, a slightly larger 
symmetrisation is needed in order to compute these ma- 
trices at the cross section level, it is interesting to remark 
that the correlation matrices at the as and cross section 
level are very similar. 



B. Results wth a simple average 

The combination using a simple average (SA) of all the 
input as values yields 



a 



SA 



n 1140+00047 

U.11'±»_Q QQ4g, 



where the quoted uncertainty is experimental only. 
These uncertainties are obtained by computing (from the 
pseudo-experiments) the variances of the positive and, re- 
spectively, negative fluctuations with respect to the nom- 
inal value. These uncertainties are rather symmetric, al- 
though the asymmetries are slightly larger when com- 
puting the values in individual (px; \y\) bins, or even the 
simple averages of the px bins in some particular \y\ bin. 



C. Results with a weighted average 

While in the previous combination method the same 
weight was assigned to all the input as values, here we 
benefit more from precise measurements, using weights 
proportional to the inverse of the squared total uncer- 
tainty. However, the gain in precision is not very big, 
because the various measurements that are used have 
rather similar precisions. Actually, all the correspond- 
ing weights are in the range between and 1, the ratio 
of the largest to the smallest weight being w 2.5. 
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FIG. 1: Examples of distributions of as{M%), as obtained from lO'' pseudo-experiments, for anti-fcf jets with R=0.6, in the 
central rapidity region (0 < \y\ < 0.3), in several pT bins. 



Computing this weighted average (WA) in the individ- 
ual \y\ intervals given in table III, we obtain 



.WA 
WA 
WA 

,WA 



y\i) 
yh) 
y\^) 
yU) 



-0.0067 

-0.0069' 

-0.0073 

-0.0077' 

-0.0075 

-0.0079' 



n 114^+00076 



0.1155: 
0.1167: 
0.1147: 



0.1147 



-0.0072 

-0.0077' 

0.0085 



v\,) = 0.1155;-- 



aWA(l2/|7) = 0.1000+S«. 

Fig. 3 shows the as{M'^) values obtained in the different 
(pt; \y\) bins, as well as the averages obtained in the cor- 
responding \y\ bins. The precision of the average in each 
\y\ bin is similar to the precision of the values obtained 
in the corresponding (input) pt bins. This is due to the 
strong correlations between (and the similar size of) the 
uncertainties in the px bins used for the average inside a 
given \y\ bin. 



Performing the combination of all the 42 (px; \y\) bins, 
one obtains 



,WA 



= 0.1151 



-1-0.0047 
-0.0047- 



(8) 



The experimental uncertainty of this average is about 
30% smaller compared to the one obtained in the most 
precise individual |y| bin (i.e. the central bin). This gain 
in precision is due to the rather small correlations of the 
systematic uncertainties, between the central and the for- 
ward region (see Fig. 2). Fig. 4 shows the global as(M|) 
average compared to the values obtained in individual 
iPT, \y\) bins and used in the combination. 

The result of this average, together with the covariance 
matrix of the uncertainties for the as values (Fig. 2(c)), 
have been inserted into Eq. 2, yielding a x^/dof 
0.54 (for 41 degrees of freedom). Using the same a 
inserted in the theoretical prediction of the cross section, 
together with the covariance matrix of the uncertainties 
for the cross section values (Fig. 2(a)), yields, by compari- 
son with the experimental cross section, a x^/dof — 0.52. 
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(b) Correlation matrix of cross section uncertainties. 
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FIG. 2: Covariance (left) and correlation (right) matrices of the cross section uncertainties (top) and as uncertainties (bottom), 
in the (px; \y\) bins used for the final as average determination. The (pt; \y\) bins are numbered as indicated in table III. 



In the individual \y\ bins, where more than one px bin is 
used in the average, the x^/dof obtained using the corre- 
sponding weighted averages and covariance matrices are 
between 0.21 and 1.39. Taking into account the bin-to- 
bin correlations is very important for the determination 
of these x^ values. Ignoring the correlations would ac- 
tually yield x^ values smaller by more than one order of 
magnitude. 

An alternative averaging procedure has been tested, 
where the weights are obtained using only the statistical 



uncertainties (which exhibit rather small correlations), 
together with the uncorrelated components of the sys- 
tematic uncertainty. However, the full information on all 
the uncertainties and correlations of the ATLAS mea- 
surement is used in pseudo-experiments, for the error 
propagation. The result obtained through this procedure 
is very similar to the previous ones (i.e. 0.1154lQ;QQ4g). 

The results (central values and uncertainties) obtained 
using the simple average and these variations of the 
weighted average are consistent, as expected for a set 
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FIG. 3: as{M'^) values (points, with asymmetric uncertainties) obtained from the ATLAS inclusive jet cross section, in all 
the available (pt; \y\) bins. The different figures correspond to increasing \y\ values, from left to right and from top to bottom. 
The green band indicates the weighted averages obtained in the corresponding \y\ bins, and it covers only the bins which are 
actually used for its determination (see text). 
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FIG. 4: Global weighted Qs(Mf ) average (green band) to- 
gether with the Qs(M|) values (with asymmetric uncertain- 
ties) obtained in all the (pt; \y\) bins used in the combination. 



of inputs with similar uncertainties. 

We have also tested the stabihty of these results when 
using non-Gaussian distributions for the experimental 
systematic uncertainties. When using log-normal distri- 
butions, the nominal values obtained from the averages 
in individual \y\ bins are very similar to the ones us- 
ing Gaussian distributions, while some small changes are 
seen in the asymmetric uncertainties. For the result of 
the combination in all the (px; \y\) bins, both the nom- 
inal value and the uncertainties are very similar to the 
ones obtained for Gaussian distributions. 



D. Minimisation of a "standard" x ^vith 
correlations 



This procedure minimises a x^ with correlation (see 
Eq. 2), using the covariance matrix in Fig. 2(c) as input, 
and yields 



0.1165: 



-0.0033 
-0.0033- 



This result has a smaller experimental uncertainty com- 
pared to the previous averaging procedure. This is not 
surprising, since the Gauss-Markov theorem guarantees 
it to have the smallest variance, among all the weighted 
averages with the sum of weights equal to unity (see sec- 
tion JVC). However, the central value obtained here is 
questionable, since the weights of this average strongly 
rely on the exact knowledge of the bin-to-bin correlations 
and are not well behaved. Indeed, a large fraction (al- 
most half) of the weights of the individual px bins in this 
average are smaller than zero. The x^/dof in this case is 
about 0.53, slightly smaller comparing to the one in the 
previous weighted average (as expected). 

Looking at the similar averages computed in individ- 
ual \y\ bins, they are in the range between 0.1100 and 
0.1243, several of them being outside the range of the 



corresponding input values. This is due to the same bad 
behaviour of the weights, some of which are negative, 
while some others are larger than one. At the same time, 
the x^/dof computed in the individual |y| bins (where 
more than one px bin is used in the average) are between 
0.20 and 0.93. 



E. Simultaneous fit of as{Mz) and the luminosity, 
using a modified x'^ 

As explained in section IV C 4, two modified x^ defi- 
nitions (at the cross section and as level, respectively) 
were introduced, treating the luminosity uncertainty as 
nuisance parameter. The main goal here is to test some- 
what better behaved definitions of the x^ with correla- 
tions. 

An approximate effective dependence of a^ on lumi- 
nosity shifts is needed in order to define the x^ ^^ the 
as level, given in Eq. 6. In order to establish this depen- 
dence, a global relative shift corresponding to one stan- 
dard uncertainty of the luminosity was performed, for all 
the cross section measurements. The first weighted av- 
erage in section V C, using the total experimental uncer- 
tainty in each (px! |y|) bin, was rederived. The induced 
variation of the weighted average yielded, within a linear 
approximation, the effective a^ (/3l) dependence (see 
Eq. 6). 

Fig. 5 shows the result of the scans of the modified 
X^ values, as a function of asiM"^) and the luminosity 
shift parameter /3l. It shows the total x^ values, together 
with their sub-components, from the comparison between 
data and theory, as well as from the luminosity shift (see 
section IV C 4). 

The partial x^ from the comparison between data and 
theory, computed at the Ofg level, exhibits a degeneracy 
between as and /3l (see Fig. 5, top left). This is ex- 
pected since, for this partial x^j a value obtained for a 
given a^ (/3l) corresponds to a combination of as and 
/3l values (cf. Eq. 6 which, for a fixed a^^ {Pl), admits an 
infinity of solutions). The degeneracy is removed when 
adding the constraint on the luminosity shifts, which fixes 
the minimum of the global x^ at /3l = 1- Therefore, the 
minimal value of the global x^ is equal to the minimal 
value of the partial x^ from the comparison between data 
and theory, and it is obtained for as{M'^) — 0.1175. One 
should recall, however, that the covariance matrix used 
here in the computation of the partial x^ from the com- 
parison between data and theory is not complete, as the 
luminosity uncertainty is treated as a separate nuisance 
parameter. This explains the difference comparing to the 
result obtained from the minimisation of the "standard" 
global x^j in the previous subsection. 

The degeneracy between as and /3lj seen for the partial 
X^ from the comparison between data and theory, com- 
puted at the as level, is absent for the partial x^ com- 
puted at the cross section level (see Fig. 5, top right). 
This is due to the slightly different dependence on as, 
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FIG. 5: x^/dof values as a function of as{Mz) and /?l (see text). The top figures siiow tiie contribution from the comparison 
between data and theory, computed at the as level (left) and at the cross section level (right). The center figure shows the 
contribution from the luminosity shifts, while the bottom figures show the total x^/dof, computed at the as level (left) and at 
the cross section level (right). 
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for the (shapes of the) theoretical calculations in the var- 
ious (pt; \y\) bins (the same effect which made Eq. 6 
to be valid only as an approximation). While the above 
mentioned degeneracy is removed, this partial x^ does 
not allow for a strong constraint of the uncertainty on 
the luminosity. After adding the constraint term and 
minimising the total x^, the uncertainty on /3l is very 
similar to the input precision of the luminosity deter- 
mination (w 3.4%). The minimum of this global x^ is 
reached by a shift of the nominal luminosity which is 
smaller than 0.5% (much smaller than the uncertainty 
on the luminosity), and an as{M^) value of 0.1160. 

While these two alternative x^ definitions still have 
some of the problems present for the "standard" global 
X^ (in particular, the negative weights obtained in the 
minimisation of the partial x^ from the comparison be- 
tween data and theory) , they show that the measurement 
of the inclusive jet cross section alone does not allow for 
an improved luminosity determination when measuring 
a,(M|). 



F. Results using a geometrical mean 

Using a geometrical average (GA) for the combination, 
we obtain: 



a 



GA 



n4q+oo"46 

U.ii^»__0.0049- 



This value and its (slightly asymmetric) uncertainties are 
very similar to the ones obtained using the simple and 
the weighted averages. The same statement is true when 
comparing the corresponding averages in the various |y| 
bins. 



G. Conclusions on the experimental results 

In the previous sections we have shown results on 
the as{M'^) determination, using several averaging pro- 
cedures. The methods using a matrix inversion in 
the x^ definition were observed to have badly behaved 
weights (negative or larger than one) sometimes yielding 
an average outside the range of the input values. Al- 
ternative averaging methods, relying on the bin-to-bin 
correlations for each component of the systematic uncer- 
tainties, have been proposed in the past. However, in our 
study, the only systematic uncertainty for which the cor- 
relations are perfectly known is the one on the ATLAS 
luminosity. We have shown that our study does not al- 
low for an improvement of the luminosity determination, 
while performing an as{M'^) measurement. 

The results obtained using various weighted aver- 
ages (with simple weights, or using the inverse of squared 
uncertainties) were shown to have (as expected) larger 
variances comparing to the ones using the minimisa- 
tion of a x^ with correlations. However, these simpler 
weighted averages have well behaved weights, which do 



not rely on the exact information on the correlations, 
while all the information on the uncertainties and their 
correlations is used in the error propagation (through 
pseudo-experiments). In addition, these results are con- 
firmed by the ones obtained using a geometrical average. 
It is for all these reasons that we choose to use the value 
obtained with the weights given by the inverse of the 
squared total uncertainties, for our nominal result. 

In Fig. 4 we have shown the nominal result for as{M'^), 
its uncertainty, and the comparison with the inputs used 
in the combination. The general agreement is good (ex- 
cept maybe for one bin from the forward region, where 
the theoretical prediction is more questionable), which is 
also confirmed by the global x^/dof (see Sec. VC). 

The N'^LO RGE running based on Runge-Kuta inte- 
gration, with 3-loop quark-fiavour matching at the top 
threshold [26, 27], has been used to evolve the as values, 
from the Z scale to a given px scale. This has been done 
for the global average, as well as for the values obtained 
in the individual (px; \y\) bins (see Fig. 6). The reduc- 
tion of the relative uncertainty on the as average, when 
performing the running from low to high px values, is 
also visible. In the same figure, we indicate the highest 
Pt value (145 GeV) used in the as{M^) determination, 
from the inclusive jet cross section, at Tevatron [28, 29]. 
Our combination procedure uses a px range going up to 
four times higher. We also use inputs from the forward 
region, up to \y\ = 4.4, while the Tevatron determination 
is restricted to \y\ < 1.6. 

The SM running of as is assumed in the various PDF 
determinations, as well as in the theoretical calculation 
of the inclusive jet cross section. A possible contribu- 
tion from "New Physics" would manifest by a deviation 
between the running of the as average and the individ- 
ual values. This would probably be more "visible" in 
the figure showing the running of the inverse of as (see 
Fig. 6) [30], although one should keep in mind that the 
same statistical information is available in all these var- 
ious ways of presenting the same result. A good agree- 
ment between our nominal fit and the 21 data points in 
the newly explored region, above 145 GeV, can be seen in 
Fig. 6 and is also shown by the corresponding partial x^j 
accounting for 8.5. The weighted average computed in 
the same region alone yields 0.1157;to.oo59, with a x^/dof 
of 0.42 (for 20 degrees of freedom). This result is com- 
patible with (but slightly less precise than) the nominal 
weighted average computed for 42 data points (Eq. 8). 

As explained above, because of the limited range in the 
as scan for the theoretical prediction, not all the (px; jyj) 
bins were used in the nominal combination. This limited 
range in the as scan induces distortions in the tails of the 
as distributions for the corresponding bins. The nomi- 
nal values are, however, less affected by this problem. 
Their combination allows for determining a systematic 
uncertainty due to the bins which where removed in the 
nominal combination, accounting for « 0.0014. 

Results with a similar precision are obtained using 
anti-fct jets with R=0.4. The central value of the cor- 
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FIG. 6: Global weighted as average (green band) evolved to the corresponding px scale (left plot) and its inverse (right plot), 
together with the values obtained in all the {px; \y\) bins used in the combination. The blue vertical line indicates the highest 
Pt value used in the Tevatron as(M|) determination, from the inclusive jet cross section. 



responding weighted average is, however, shifted down- 
wards with respect to our nominal result (from anti-fct 
jets with R=0.6), by about 0.0060. This is expected, 
due to a systematic difference seen between the compar- 
isons of the two (strongly correlated) measurements with 
the corresponding theoretical predictions (see section II) . 
The study of the x^ scan (see section V E) applied to anti- 
kt jets with R=0.4, indicates an even smaller "preferred" 
shift of the luminosity value, comparing to the one ob- 
served for anti-fct jets with R=0.6. It should be pointed 
out that a shift in the nominal value of the luminosity (or 
of any other nuisance parameter which is strongly corre- 
lated between the two measurements) could not explain 
the observed difference anyway, as these shifts should be 
identical for anti-fct jets with R=0.4 and R=0.6. The 
difference between the ag determinations with the two 
jet sizes is treated as systematic uncertainty in our final 
result. This is actually our largest uncertainty and an im- 
provement in the understanding of these differences is de- 
sirable. This improved understanding could be achieved, 
for example, through the measurement of the ratio be- 
tween the spectra of anti-A;* jets with R=0.4 and R=0.6, 
and the comparison with the corresponding theoretical 
prediction (as given in [18]). The main challenge of this 
measurement is the understanding of the (strong) corre- 
lations of the systematic uncertainties of the two spectra. 



{(i,l), (2,1), (l,i), (l,2),(i,i), (2,2)} in each 
[pt, \y\) bin. The nominal value of the measured 
cross section was mapped to the ag (/j,f, Mfl) for each 
scale choice. The asymmetric envelope of the values 
as {fJ-F^ (J-r) around the one with the default scale choice 
is taken as the scale uncertainty on as- Propagating this 
uncertainty through our nominal averaging procedure 
results in a scale uncertainty on the as{M'^) value of 
-h0.0044 -0.0011. 

For the evaluation of the uncertainty on as due to the 
choice of the PDF set, the theory prediction for the inclu- 
sive jet cross section was recalculated using MSTW 2008, 
NNPDF 2.1 (100) [8, 9] and HERAPDF 1.5 [10] PDF 
sets. For each alternative PDF set the a^ ^^'"^* was de- 
termined from the measured cross section and the recal- 
culated theoretical one, using the procedure described in 
Section IV. The estimated value of the PDF choice un- 
certainty is -1-0.0022 — 0.0015, which is the maximal dif- 
ference between the af ^^'''^* and the central CTIO PDF 
one. 

For the CTIO PDF the A'^p = 26 parameters (eigenvec- 
tors) were determined. Therefore CTIO provides 2 x Np 
PDF uncertainty sets for positive and negative varia- 
tions of each eigenvector in the Ap-dimensional param- 
eter space [5]. The asymmetric PDF uncertainty on as 
for the CTIO PDF is given by [31]: 



H. Theoretical uncertainties 

The main uncertainties on the NLO QCD prediction 
come from the choice of the PDF set, the uncertainties on 
the nominal set of PDFs and the choice of factorisation 
and the renormalisation scales. 

In order to propagate the scale uncertainty on 
the cross section to the uncertainty on as, the 
theoretical cross section was recalculated with 
six different scale choices (- 



Aat - 



Ac 



N„ 



\ X! t™^^ (("s )i - (as)o, {as )i - (as)o, O)] ' 



N„ 



\ X! t™'^^ ((as)o - {as)t, (as)o - ("s), , O)] , 



where {as)o and {af)i are calculated using the central 
PDF set and ± variation of all PDF set eigenvectors, re- 
spectively. The PDF uncertainty on the average as{M'^) 
accounts for ±0.0010. 
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Finally, for the evaluation of the uncertainty on as due 
to the modelling of non-perturbative effects, the nomi- 
nal value of the measured cross section was mapped to 
the a^*^*"'" for each MC tune used to estimate the un- 
certainty of the non-perturbative corrections on the in- 
clusive jet cross section. Again, the asymmetric enve- 
lope of the values 0;^^*^*""= around the one with the de- 
fault tune is taken as the uncertainty on as{M'^) due 
to the modelling of non-perturbative effects, yielding 
-1-0.009 -0.0034. 

Given the fact that the theoretical systematic uncer- 
tainties from the scale choice, PDF eigenvectors, the 
choice of the PDF set and modelling of non-perturbative 
effects, are (to a good approximation) independent, the 
total theoretical uncertainty is equal to their quadratic 
sum. 



where the uncertainties are statistical and experimental 
systematic (propagated from the ATLAS data), due to 
the limited px range, due to the differences between the 
results obtained with the R=:0.4 and R=0.6 jet sizes, due 
to the renormalisation and factorisation scale choice in 
the theoretical calculation, due to the choice of the PDF 
set, propagated from the PDF eigenvectors, and due to 
the non-perturbative corrections, respectively. Our value 
is in good agreement with the latest (preliminary) update 
of the as(M|) world average (0.1183 ± 0.0010) [32], as 
well as with the latest result from a hadron-hadron col- 
lider (O.lieitaools) [28, 29]. Although our result is less 
precise, it includes for the first time the measurements of 
the inclusive jet cross section up to 600 GeV. The run- 
ning of as has also been tested, in the pt range between 
45 and 600 GeV, and no evidence of a deviation from the 
QCD prediction has been observed. 



VI. CONCLUSIONS 

We have performed a determination of the strong cou- 
pling constant at the Z scale, using the ATLAS inclusive 
jet cross section data. Our final result accounts for 



a.(M|) 



0.1151 ± 0.0001 (stat.) ± 0.0047 (exp. syst.) 
±0.0014 (pT range) ± 0.0060 (jet size) 
t°oTit (scale)lHSf5 (PDF choice) 
±0.0010 (PDF eig.)l°;°°°^ (NP corrections). 
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